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RAREFACTION  SHOCK  POSSIBILITY 
IN  A  VAN  DER  WAALS -MAXWELL  FLUID 


ABSTRACT 


The  rarefaction  shock  is  found  to  be  impossible  practically  in 
an  equilibrium  Van  aer  Waals-Maxwell  fluid,  in  and  across  all  phase 
regions.  The  existence  of  this,  type  shock  depends  explicitly  on  the 
const ant -volume  specific  heat,  as  well  as  the  pressure  equation  of 
state.  The  form  of  specific  heat  assumed  here  includes  contributions 
from  translation,  rotation,  and  vibration  energies  of  the  molecule; 
Einstein  functions  are  used  to  represent  the  latter.  In  particular, 
the  vaporizing  expansion  wave  in  the  saturated  liquid  phase  cannot  be 
discontinuous.  Metastabie  states  of  super saturation  are  not  considered. 
Despite  the  results  for  this  theoretical  model,  the  necessary  conditions 
for  the  rarefaction  shock  are  found  to  be  satisfied,  in  principle,  for 
a  small  region  near  the  critical  point  of  a  real  two-phase  fluid 
(steam) . 
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I.  INTRODUCTION 


In  most  fluids,  including  the  ideal  gas,  an  initial  compressive 
disturbance  steepens  with  time  into  a  compressive  shock  discontinuity, 
while  an  expansive  disturbance  flattens  out,  never  producing  a  dis¬ 
continuity.  This  type  of  response  depends  on  the  equation  of  state  for 
the  fluid.  In  principle,  the  possibility  of  contrary  responses  to 
initial  disturbances  exists  in  fluids  obeying  certain  state  equations. 
For  example,  Paxton^  discusses  the  rarefaction  shock  possibility  in  a 
plasma  carrying  a  large  electrical  current.  And  Zeldovich  and  Raizer2 
cite  a  possible  rarefaction  discontinuity  in  the  gaseous  phase  of  a 
Van  der  Waals-Maxwell  fluid  near  its  critical  point;  their  treatment  of 
this  problem  is  a  reproduction  of  one  given  earlier  by  Zeldovich3  in 

19^7. 

A  detailed  analysis  of  the  shock  discontinuity  for  arbitrary  state 
equations  was  made  in  19^2  by  Bethe4  in  a  report  of  initially  limited 
distribution.  And  it  is  shown  in  [4]  that  the  probability  of  satis¬ 
fying  a  necessary  condition  for  the  non-MHD  rarefaction  shock  in  the 
gaseous  phase  of  a  substance  is  practically  non-existent,  including 
the  specific  fluid  treated  in  [2]  and  [3].  The  equilibrium  coexistence 
phase  of  both  liquid  and  vapor  was  not  considered. 

In  this  paper,  the  two-phase  region  of  the  Van  der  Waals-Maxwell 
fluid  is  examined.  Motivation  for  this  investigation  is  provided  by 
observations  on  exploding  wires.  The  early  electrical  resistance  of 
many  exploding  wires  fits  a  model  which  assumes  a  vaporizing  expansion 
wave  proceeding  radially  inward  from  the  cylindrical  surface,  shrinking 
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the  conducting  cross  section  of  the  intact  molten  interior  ’ 


The 


existence  of  such  a  wave  can  be  predicted  from  a  two-phase  fluid  model; 
and  the  Van  der  Waals-Maxwell  fluid  yields  wave  velocities  which  are  in 
fair  agreement  with  experiment7 .  It  is  of  interest  to  know  whether 
this  model  will  predict  a  steepening  of  this  expansion  wave  into  a 
finite  rarefaction  shock.  Although  the  discontinuous  part  of  the 


*References  are  found  on  page  29. 
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expansion  is  not  an  absolute  necessity  for  the  experimental  interpre¬ 
tation,  it  would  assist  in  offering  a  simple  and  conceptually  plausible 
explanation  for  the  resistance  behavior. 

This  treatment  is  limited  to  thermodynamic  equilibrium,  thereby 


mrnl  i  +-  c  one  en  /->Vi  o  o  c,nT>av»ViooH — i  ran-  t~\  m 

CAL.xuuin5  mciaouaui-c  LU  liU.  _L  L  lL-iio  ou  11  x  nv- a,  U  .Llig,  Ui 
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one  of  the  phases.  Although  ever  mounting  evidence  shows  the  Van  der 
Waals -Maxwell  model  is  not  an  accurate  description  in  the  critical 
region  of  a  fluid,  it  is  still  valuable  for  its  simplicity,  and  as  a 
standard  of  comparison  for  recent  empirical  proposals  of  state 
functions  there.  Recent  computations  of  its  coexistence  properties 
have  been  made  by  Barieau8 ,  and  we  adopt  his  notation  for  the  scaled 


We  find  that  the  rarefaction  shock  possibility  is  even  less 
likely  in  the  two-phase  region  of  this  fluid  than  in  the  single  phase, 
where  it  has  already  been  shown  to  be  practically  impossible4.  Then, 
despite  the  results  with  this  model,  we  show  by  direct  use  of  steam 
tables,  as  suggested  by  Novikov9,  that  the  rarefaction  shock  possibil¬ 
ity  cannot  be  dismissed  over  a  small  part  of  the  wet  vapor  state  of 
steam  near  critical. 


II.  DISCONTINUITY  CONDITION 

Let  P,  V,  T,  and  S  be  pressure,  specific  volume,  temperature,  and 
specific  entropy  respectively.  The  change  of  entropy  across  a  dis¬ 
continuous  transition  of  a  system  initially  at  the  state  (P0 ,  V0 ,  T0 , 
S0  )  to  the  state  (P,  V,  T,  S)  is  given  by 

s  -  So  =  -  (V  -  VO)3{(0)  }  /(12To)  +  •••.  (1) 
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This  expression  follows  from  several  properties  concerning  the 
mechanical  discontinuity  relations  which  are  proved  by  Courant  and 
Friedrichs10.  These  properties  are:  a)  the  pressure  function  P  =  G(v) 
of  the  transition  Hugoniot  curve  makes  a  second  order  contact  with  the 
pressure  function  P  =  g(v)  of  the  adiabatic  isentrope  at  the  conditions 
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of  the  initial  state;  and  b)  the  entropy  change  along  the  Hugoniot 
curve  is  third  order  in  volume  change  from  the  initial  state  with 

(§)„  =  (|p)?  *  °>  an®  (§f\  =  -  (3 The  first  condition 

implies  1  (rnl-)  ]■  =  f— §-)  ;  making  a  series  expansion  of  (S-S0  )  in 
(AoV“  1  s  .1 0  ^ CIV''  / o  ~ 

terms  of  (V-V0 )  then  results  in  the  first  non-zero  term  shown  in  Eq. 
(l)«  Not  all  gas  dynamic  transitions  possess  the  foregoing  simple 
properties;  a  notable  exception  is  the  detonation  wave.  For  our 
purpose,  we  have  assumed  that  the  fluid  on  both  sides  of  the  disconti¬ 
nuity  obeys  a  unique  state  equation,  so  that  the  transition  is  neither 
exothermic  nor  endothermic;  a  discussion  of  the  Hugoniot  curves  for 
these  latter  transitions  is  given  by  Hayes11 . 

The  form  of  Eq.  (l)  shows  that  for  the  usual  fluid  where  (=^?r) 

“  >  •  s 

positive,  the  transition  which  increases  the  entropy  forces  V  <  V0 ,  as 

/a2  p\ 

in  a  shock  compression.  Were  negative,  however,  then  the 

condition  for  increased  entropy  would  be  V  >  V0 ,  corresponding  to  the 
rarefaction  discontinuity.  The  possibility  of  the  rarefaction  shock 


process  is  seen  to  depend  on 


(|y?)  being  negative. 


It  is  this  quanti¬ 


ty  which  we  investigate  for  the  Van  der  Waals-Maxwell  fluid.  The  state 
functions  P,  E,  and  S  (where  E  is  the  internal  energy)  are  expressed  in 
terms  of  the  independent  variables  V  and  T.  We  will  use  the 
differential  operator 


(*-)  =  (p)  (*)  +(p) 

\dV/8  V9V/S  \9T/V  VdV/T 


(2) 


along  with 


TdS  =  dE  +  PdV 


(3) 


and  the  identity 


Both  dE  and  dS  are  assumed  exact  differentials  of  V  and  T  within  given 


V,  T,  regions.  The  constant  volume  specific  heat  is  Cv  =  (I).  ,  and  i : 
assumed  never  negative.  From  Eq.  (3)  one  finds  =  -  P  directly, 

and  application  of  operator  (2)  to  E,  and  using  (4)  results  in 
(H)  =  "  <T/Cv  )  (|?)  ’  a§ain  using  Eq.  (f)  on  P  results  ii 


in 


(as).  -  -  <*/c. )  (n):  ♦  (i), . 


(5) 


With  w  the  adiabatic  (and  in  this  instance,  isentropic)  sound  speed  of 

the  substance  defined  by  w2  =  -  V2  (^)  ,  the  right  side  of  Eq.  (5)  is 

seen  to  be  -  w2 /V2  .  The  sound  speed  is  expected,  to  be  real  and  finite 
for  a  well-behaved  fluid;  sin^e  T  and  Cv  are  always  positive,  then  w 

dP\ 


will  certainly  be  real  if  \^yj  s  0.  This  latter  condition  holds  for 
most  substances,  including  the  Van  der  Waals-Maxwell  fluid,  which  we 
now  describe  more  completely. 


III.  FLUID  MODEL 

The  fluid  may  be  either  a  liquid,  vapor,  or  an  equlibrium  mixture 
of  the  two  phases.  For  either  of  the  single  phases,  the  pressure  is 
given  by  the  Van  der  Waals  function, 

PW(V,T)  =  [rt/ (V-b )}  -  a/v2,  (T  >  0,  V  >  0)  (6) 


where  R  is  the  gas  constant  for  unit  mass  and  a  and  b  are  constants  typi¬ 
fying  the  substance.  By  choosing  for  unit  mass  the  molar  mass  of  the 
substance,  one  replaces  R  with  R,  the  universal  gas  constant.  The  (V,T) 
region  where  Eq.  (6)  applies  is  T  >  0,  V  >  0  except  for  the  region  of 
phase  coexistence .  The  coexistence  region  is  defined  by  T  ^  Tc  (where 
subscript  c  always  denotes  critical  conditions)  and  V3(t)  ^  V  s  v  (t), 


1? 


with  the  loci  V3(t)  and  Vi (T)  specified  later.  Within  the  region  where 


+  hn+  r 


Eq.  (6)  applies,  a  differentiation  by  T  of  Eq.  (4)  shows  „v 

independent  of  V.  Outside  coexistence,  the  fluid  consists  of  a 
single  phase,  either  liquid  or  vapor.  For  convenience  we  label  the 
regions  T  ^  Tc ,  b  ^  V  ^  V3 (T)  as  the  liquid  phase,  and  assert  the  vapor 
phase  occurs  in  both  T  ^  Tc ,  V  ^  b,  and  in  T  s  Tc ,  Vi (l)  s  V. 

The  pressure  in  the  equilibrium  coexistence  region  is  not  PW(V,T), 
but  is  instead  the  vapor  pressure  function,  Pa(t),  independent  of  V. 
This  function  is  not  defined  explicitly;  it  must  satisfy  two  conditions 

a)  the  pressure  across  the  coexistence  boundary  must  be  continuous,  and 

b)  the  chemical  potential,  or  Gibbs  free  energy  function  E  -  TS  +  PV, 
is  constant  for  fixed  T  within  coexistence,  and  is  continuous  across  th 
boundaries.  Both  these  conditions  are  satisfied  by  the  Maxwell  rule: 
for  any  T  ^  Tc  , 


Vi(T) 

■* 

Pw(V,T)dV  =  Pa(T)  (Vi(T)  -  V3  ( T ) }  ,  (7) 

i 

V3(t) 

with  the  subsidiary  conditions  P«  (Vi  ,T)  ’=  PA  (T)  =  Pw(V3,T).  Eq.  (7) 
and  its  subsidiary  conditions  fix  PA (T)  and  the  V2 (t)  and  Vi  (t)  loci  of 
the  coexistence  boundary.  This  Maxwell  rule  implicitly  imposes  the 
restriction  for  fixed  T  ^  Tc  that  Cv (t) (liquid)  =  Cv(T)(vapor)  =  cj (T) . 
Although  there  is  another  suggested  rule  similar  to  Eq.  (7)  which 
relaxes  this  condition  on  the  specific  heats12 ,  it  is  less  convenient 
because  it  requires  additional  information  about  the  difference  of 
these  specific  heats.  For  simplicity,  we  confine  present  attention  to 
the  Van  der  Waals -Maxwell  fluid  with  the  coexisting  phase  region 
governed  by  Eq.  (7). 

The  isotherm  segments  of  PA (t)  are  isobaric  and  extend  from  V3(t), 
the  saturated  liquid  locus,  to  Vj(t),  the  saturated  vapor  line.  By 
replacing  those  portions  of  the  mathematically  defined  but  thermo¬ 
dynamically  unacceptable  isotherm  segments  of  Eq.  (6),  with  the  accept¬ 
able  isotherm  segments  of  PA (t)  in  the  coexistence  state,  we  have 
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/?>p\ 


removed  any  possibility  that  might  be  positive  for  this  fluid. 

This  fact  insures  that  the  right  hand  side  of  Eq.  (5)  is  negative  in 
all  regions  of  interest. 


IV. 

We  review  briefly  the  conditions  found  in  the  single  phase  of  thi 
fluid.  Allowing  Cv  to  be  a  function  of  temperature  alone,  and 
successively  applying  operator  (2)  twice  to  Eq.  (6),  one  finds 


( —  ■  fid  ^  _  f  dt /  (\t  i.  ^  1  v/ir  ml] 


) 


with  the  dimensionless  quantities 


0 (T)  =  (2  +  R/C° )(1  +  R/C° ) 


T(R/c;  f  fpv 


(9) 


Y(V,T)  =  6a(V-b)3 / (RTV4 ) . 


(10) 


We  see  two  ways  the  right  hand  side  of  Eq.  (8)  can  be  negative. 

First,  if  0(T)  <  0,  the  right  side  is  always  negative  because 
Y(V,T)  is  >  0,  and  V  £  b.  For  0(TJ  <  0,  the  rate  of  changes  of  C° (t) 
with  increasing  T  must  be  large.  If  0(T)  <  0,  then  even  the  ideal  gas 
might  sustain  a  rarefaction  shock.  This  fact  may  be  seen  by  observing 
that  the  Van  der  Waals  pressure  function  reduces  to  that  of  the  ideal 
gas  when  both  constants  a  and  b  of  Eq.  (6)  are  identically  zero;  the 
right  side  of  Eq.  (8)  then  becomes  0(t)RT/V3,  with  sign  depending  on 
0(t).  Thus,  negative  0(t)  permits  both  the  ideal  gas  and  the  single 
phases  of  the  Van  der  Waals  fluids  to  meet  the  necessary  condition  for 
the  rarefaction  shock.  The  probability  of  negative  0(T)  will  be 
discussed  in  Section  VII. 

But  there  is  a  second  way  that  the  right  side  of  Eq.  (8)  might  be 


negative.  This  condition  occurs  if 


Y(V,T)  >  0(t) 


(11) 


for  a  particular  (V,T)  domain.  The  constants  a  and  b  can  be  expressed 
in  terms  of  the  critical  values  of  P,  V,and  T,  as  b  =  Vc/3, 

~  _  fr\/Q\-n  m  it  - —  1  ^4-  Tb  it  _  ( o  /Q  "Norn  nru^  -p,^  ^4-  A  V  At  rp  \ 

ct  —  \i^/^^rvxcvC  ctllu-  '--lie  rexaulun  rc  Vc  —  xnc  x  luiu  G  iun  j.  \,  v  ^  x  / 

has  a  maximum  value  with  respect  to  V  at  V  =  4b  =  (4/3  )VC  ;  meanwhile 
the  permissible  T  values  are  bounded  below  by  T(min),  corresponding  to 
values  along  "the  coexislence  "boundary $  since  Ec[«  (8)  applies  only  "to 
the  single  phase  region.  At  V  =  (4/3 )VC ,  T(min)  =  .9834  Tc ;  if  Cv  is 

constant  with  T,  so  =  0,  then  inequality  (ll)  can  be  satisfied 

for  the  dimensionless  quantity  (C°/R)  ^  17. 5-  Very  large  values  of 
-  o  cKCu  /r'' 

(Cv/r)  and/or  -  are  needed  to  fulfill  inequality  (11),  even  when 

the  specific  volume  lies  in  a  small  neighborhood  slightly  larger  than 

o  , 

critical  volume.  Again  we  defer  discussion  of  Cv/R  to  a  later  section. 


V.  COEXISTENCE  STATE 
A.  Justification  of  Model 

This  section  deals  only  with  the  (V,T)  region  of  coexisting  phases 
of  the  fluid.  However,  some  justification  and  caution  is  needed  in 
applying  continuum  fluid  mechanical  concepts  to  this  regime.  The 
vapor  pressure  function,  PA (t) ,  used  in  the  equilibrium  two-phase  region 
applies  to  a  variety  of  combinations  of  coexisting  liquid  and  vapor 
phases,  requiring  only  that  the  phases  have  equal  temperature  and 
pressure.  For  example,  the  entire  liquid  phase  could  be  located  at  one 
end  of  a  container,  with  the  vapor  phase  at  the  other;  in  this  instance, 
a  unique  "sound"  speed  for  the  system  would  appear  unjustified.  The 
two-phase  fluid  mixture  could  approximate  a  continuous  medium,  pro- 

in  rli  nrr  oo  nVi  ttVi  o  o  o  to  rl  *i  ono  ■poor]  t.n  4-  V  i  xi  4-  V*  r\-f  Viq-p  o/t  4-  Vi  4-  4-  Vi  w?  v4-  4  « 

v  xuj.  ^  u.^n  j-v  li-cxu  -i_  o  iu.  j.  x  o  v-i.  w  _l  uuxii  wu'C  u  ouci  o  '  onai<  O  uc  iiiJ-A  o  * *  t  C?  _L11 

a  small  volume  appears  homogenous  on  the  length  scale  of  that  volume. 
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.ddition,  the  dimensions  of  this  volume  are  small  enough  that 


the  mechanical  relations  may  be  reasonably  approximated  by  differential 
equations,  then  continuum  fluid  mechanics  can  be  used.  We  assume  that 
these  conditions  apply  in  our  further  treatment  of  the  coexistence  state, 
In  this  connection,  Cowperthwaite  and  Ahrens13  have  recently  given  a 
discussion  of  the  thermodynamics  of  the  adiabatic  change  of  a  system 
consisting  of  two  phases  at  the  same  pressure.  They  find  several  ways 
such  a  change  is  also  isentropic;  one  of  these  requires  reversible  heat 
exchange  between  the  phases,  with  the  mixture  composition  changing  so 
that  the  phases  remain  in  chemical  equilibrium.  Using  the  vapor 
pressure  function,  Pa(t),  defined  by  Sq.  (?)  as  the  pressure  state 
equation  for  the  coexistence  region,  insures  equilibrium  of  both 
chemical  potential  and  pressure;  thus,  infinitesimal  adiabatic  changes 
in  this  two-phase  system  will  also  be  isentropic.  The  question  of 
whether  the  finite  adiabatic  change  of  this  type  remains  isentropic 
requires  further  investigation,  and  it  is  this  problem  that  we  are 
considering, 

B.  Specific  Heat 

By  using  PA  (rl)  for  P  in  Eq.  (4),  integrating  with  respect  to  V 
from  V3(t)  to  V(t),  where  (v3  (T)  ^  V  s  Vi(t)},  and  differentiating  the 
result  by  T,  one  obtains  the  specific  heat  at  constant  volume  for  the 
coexistence  region.  It  is  CvA(t),  given  by 


Cva(V,T)  =  Cv(T)  +  A(T)  +  (V-V3  (t)}  J(t), 


\lc-  I 


where  A(t)  =  V3'(PA-TPA'  +  a/V33 )  and  J(T)  =  TPA " ,  where  the  prime 
hereafter  signifies  ordinary  temperature  derivative,  and  Cv (t)  is  the 
constant  volume  specific  heat  outside  coexistence.  As  expected, 
Cva(v,T)  is  a  linear  function  of  V.  A  further  useful  relation  can  he 
obtained  from  Eq.  (7)?  viz., 


TPa  PA  =  a/(v1V3). 


(13) 
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Eq.  (13)  results  from  differentiating  Eq.  (7)  by  T,  applying  the  sub¬ 
sidiary  conditions,  multiplying  the  result  by  T  and  subtracting  Eq,  (7) 
from  both  sides.  Computed  values  from  Eq.  (7)>  as  well  as  limiting 
values  as  T  -1  Tc  *  show  that  PA  ' ,  PA  ",  and  V3  '  are  all  positive  for 
t  s  Ic ,  It  follows  that  j(T)  is  always  positive  in  coexistence;  since 
Vx  ^  V3 ,  using  Eq.  (13)  in  the  definition  of  A(t)  shows  that  A(t)  is 
never  negative  here.  Therefore,  CvA (V,T)  is  positive;  since 
Vi  >  V  ^  va  for  any  fixed  T,  one  finds  the  inequality 


0  <  CvA(V3,T)  ^  Cva(V,T)  ^  CvA(Vi,T).  (14) 

C .  Coexistence  Isentropes 

Since  PA (T)  is  independent  of  V;  using  PA (t)  for  P  in  Eq.  (5) 
gives  =  ~  (t/CvA  )(P*  ')2  *  Using  operator  (2)  on  this  latter 

expression,  remembering  that  CVA  depends  on  both  T  and  V,  results  in 


(!^-)s  =  (t/Cv  A3  ) (PA  '  )3  L(V,T), 

where 

L(V,T)  =  {l+3J(T)/PA'3  Cva  -  T  (^-)v  • 


From  Eq.  (12)  one  finds 


where 


(15) 

(16) 


h7^ 
\  i  ) 


(18) 


■in  PTyDPnHiv  +VIQ+  TC  nncji  4-n  irp  ■frjT1  rP<T.  W"1  "t*  V) 

V  ^  V3  (t)  in  coexistence,  the  following  inequalities  come  from  Eq.  (ly) 


Noting  that  PA '  and  J(T)  are  positive  quantities,  one  ran  make  use  of 
relations  (l4),  (19),  (l6)  and  (15)  to  obtain  the  further  inequality 


L(v'i  ,T)  ^  L(V,T)  s  L(V3,T), 


(  f~\r\  \ 


and  consequently 


{(IS1-)  (Vl  ,T)1  <  {(IS1)  (V,T)}  *  {(P1-)  (V,,T)} 


(?l  'i 


Ineauality  (21 )  shows  that  the  size  of  the  second  volume  derivatives  of 
the  isentropes  in  coexistence  is  bounded  between  the  value  on  the 
saturated  liquid  and  vapor  lines,  the  smaller  being  on  the  latter  locus. 
Thus,  examination  for  the  rarefaction  shock  possibility  can  be  made  at 
the  saturated  vapor  locus,  with  the  assurance  that  if  the  necessary 
condition  be  met  anywhere  in  coexistence,  it  must  also  be  met  on  this 
locus.  We  recall  from  Section  IV  that  the  most  likely  neighborhood  for 
the  occurence  of  this  condition  in  the  single  phase  was  also  adjacent 
to  the  saturated  vapor  line,  i.e.,  for  V  «  (4/3)  Vc  and  T  =  T(min). 

By  Eq.  (15),  the  sign  of  (|^ I4-)  is  seen  to  be  that  of  L(V,T); 

setting  V  =  Vi (T)  in  Eq.  (l6)  gives  L(Vi ,T)  on  the  saturated  vapor 
locus.  Defining  Z(T)  =  1  +  3J /P*  ' ,  and  arranging  the  terms  to  show 
explicitly  the  dependence  on  the  specific  heat  of  the  single  phase,  one 
obtains 


with 


L(Vi  ,T)  =  ZC°  -  T(C°)'  +  K(T), 


(22) 


k(t)  =  zm(t)  -  tn(t), 


W) 


M(T)  =  Cva(Vi,T) 


C°  =  A  +  [Vi-V33  j(t)  , 


(24) 


-  r> 


lO 


and 


N(T)  =  |~|(CV  A~P°,).^ 


(Vi,T)  =  A  -  JV3  '  +  (V1-V3)  J' 


V^) 


All  the  terms  making  up  the  function  K,  M,  and  N  are  explicit  functions 
of  T,  independent  of  C°  ,  and  can  be  computed  numerically  from  Eq.  (7). 


Hit:  liuiiicx 


ical  methods  are  satisfactory,  except  when  T  approaches  Tc ; 


but  here  the  limits  of  the  functions  can  be  found  from  a  series 
expansion  around  Tc8 .  The  computations  are  conveniently  made  in  terms 
of  the  reduced  variables  formed  by  scaling  P,  V,  and  T  by  their 
critical  values. 

The  limits  at  Tc  of  most  of  the  functions  making  up  Z(T)  and  K(t) 
are  given  in  Ref.-  [8 ] 5  the  few  others  needed  are  shown  in  our  appendix. 
The  forms  of  Z(T)  and  k(t)/R  are  displayed  in  Figs.  1  and  2;  they  are 
both  always  positive,  monotone  decreasing  with  T,  and  have  their 
smallest  values  of  8.200  and  31-86  R,  respectively,  at  T  =  Tc .  With 
Z  and  K  always  positive  in  coexistence,  Eq.  (22)  shows  that  if  C$  is 
constant  with  respect  to  T,  t (C°  )  =  o] ,  there  is  no  chance  of  a  rare¬ 

faction  shock  in  the  two -phase  region.  This  fact  contrasts  with  the 
the rmo dynamic  behavior  of  this  fluid  in  the  adjacent  vapor  phase,  where 

fv*r\ 

a  constant  but  large  value  of  Cy  permits  negative.  The 

only  possibility  of  negative  L(V,T)  rests  with  T(Cy ) '  being  larger  than 
ZC°  +  K(t),  or  in  the  form  of  an  inequality, 


T(C$/R)  ^ [K(T)/R]  +  Z(T)  (C°/r)&  8.2  (Cy/R)  +  31-86  . 


(26) 


)  across  the  coexistence 


D .  Coexistence  Boundary 

For  completeness,  the  change  of  \^j 

/  “N  T— .  \ 

boundary  should  be  investigated.  Thus,  a  comparison  of  (^k1-)  for  the 

WV  '  e 

/g  p  \ 

single  phase  isentrope  with  V5v7.  of  coexistence,  each  evaluated  at 
the  same  boundary  point  is  needed  to  show  the  sign  of  any  discontinuity 


3  Q 
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Figure  2.  The  dimensionless  function  K(t)/R  (See  Eqs.  (23-25)} 
vs  T/Tc . 
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in  (^)  across  this  boundary.  This  problem  on  the  phase  boundary  for 

the  liquid-vapor  transition  has  previously  been  treated  in  Ref.  [4], 
for  an  arbitrary  equation  of  state;  the  discontinuity  was  always  found 
to  be  positive.  The  Van  der  Waals -Maxwell  fluid  is  no  exception  to 

this  rule.  We  recall  from  Section  II  that  is  everywhere  negative 

for  this  fluid.  Omitting  considerable  algebra  we  exhibit  the  result 

for  A  =  -  (gy*")  5  where  ail  quantities  are  evaluated  at  a  point 

on  the  phase  boundary,  either  saturated  liquid  or  vapor;  the  product  of 
A  and  Cv  A (Vj  ,TJ  is 

4C,.(V,,T)  =  -  (l/c!)  (|Sl)  {c;  +  tv3'(|Sl)  f  -  fj  (|^)  ,  (27) 


with 


fj  =  TPa  "  {Vj  (T)  -  V3(t)} 


(28) 


and  subscript  J  taking  the  value  1  or  3>  depending  on  whether  the 
boundary  point  lies  on  the  saturated  vapor  or  liquid  locus,  respec- 

ively.  The  quantity  f,  is  either  positive  or  zero,  while  is 

always  negative,  so  the  contribution  from  the  second  term  of  Eq.  (27) 
is  either  positive  or  zero  because  of  its  negative  sign.  Meanwhile 


/^xtj  \ 
/  UJTf  * 


'§V  ^  is  never  positive  on  and  outside  the  coexistence  boundary, 
while  Cv A  is  always  positive,  so  A  is  never  negative.  Hence,  no  rare¬ 


faction  discontinuity  can  occur  because  of  a  negative  jump  in  (v 
across  the  coexistence  boundary. 


ria 


dV/  < 


VI.  SPECIFIC  HEAT  FUNCTION 

For  this  fluid,  a  necessary  condition  for  the  rarefaction  shock 
depends  explicitly  on  the  constant  volume  specific  heat.  We  explore 
briefly  the  expected  contribution  from  a  representative  specific  heat 
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«  n 

function.  Choosing  the  same  units  for  Cv  and  R,  the  ratio  Cy /R  is 
dimensionless.  To  good  approximation,  Cj /R  can  be  expressed  in  terms 
of  translational,  rotational  and  vibrational  energy  contributions, 
respectively,  as14, 

(C?/R)  =  (C?/R)(tr.)  +  (C°v/R)(rot.)  +  (C°  /r)  (vib . ) .  (29) 

T.T~ 

w  6  cx  tune 

(C°  /R)  (vib . )  gt  (xt),  (30) 

1 

where  all  sums  range  from  1  to  A,  where  A  is  the  number  of  vibrational 
modes,  and  g,  (x,  )  is  an  Einstein  function  given  by 


gt  (*i  )  =  (xj  /sinh  Xj  )2  , 


(31) 


with  Xj  =  9j/(2T),  and  0,  a  constant  for  the  i^h  mode.  Eq.  (30)  is 
usually  applied  to  a  polyatomic  gas,  where  the  first  two  terms  on  the 
right  sum  to  a  constant  value  B/2  ranging  between  3/2  and  3,  according 
to  the  molecular  structure,  and  A  =  3N  -  B  where  N  is  the  number  of 
atoms  per  molecule.  If  N  >  2,  B  is  always  6,  a  condition  we  henceforth 
assume.  If  optical  modes  (electronic  vibration)  are  also  included,  A 
becomes  correspondingly  larger,  as  N  must  now  include  the  electrons; 
however,  the  9,  for  the  optical  modes  are  usually  so  high  that  they  add 
little  to  the  specific  heat  at  temperatures  less  than  5000°  K1 4 . 

From  the  forms  of  Eqs.  (29)  and  (30),  it  follows  that 


ncy/R 


A 

) '  =  Y,  kt  ) 


(32) 


with 


\ 

1 


-  2  r„ 


A. 


-  s  inh  Xj  ]/  (  s  inn  xt  ) 


\3 


°3 


This  function  kj (xt )  has  a  maximum  value  of  .6757  at  xt  =  1.744;  while 
g.  (xj  )  has  its  largest  value  of  1  at  x,  =  0.  Therefore, 

T(C?,  /r)'  £  .676  A,  (34) 


where  equality  is  approached  only  if  all  0  have  the  same  value,  and 
k,  is  at  its  maximum;  and 


(35) 


where  the  equality  is  approached  only  if  every  8  4  T,  so  each  x1  is 

nearly  zero.  With  the  assumption  B  =  6,  use  of  Eas,  (29)  and  (35)  give 

3  *  C°/R  ^  3  +  A  .  (36) 


VII.  RAREFACTION  SHOCK  POSSIBILITY 


With  the  preceding  model  for  the  specific  heat,  one  may  now  assess 


the  chance  of  satisfying  the  necessary  condition 


<  0  for  the 


rarefaction  shock  in  this  fluid. 


A .  Ideal  Gas 

One  seeks  conditions  where  0(t)  might  be  negative.  Making  use  of 
the  strongest  form  of  the  inequalities  Eqs.  (34)-(36)  in  the  definition 
of  0(T),  yields  the  very  strong  inequality 


0(T)  >  2  -  .676  A/27  . 


(37) 


One  concludes  that  0  could  not  possibly  be  negative  for  A  less  than  80 
modes:  this  condition  appears  to  exceed  reality  sufficiently  as  to  be 


practically  impossible.  Actually,  a  weaker  form  of  inequa 


lity  indicate 


A  must  be  very  much  larger,  as  pointed  out  in  Ref.  [4].  For  example, 
assuming  all  0  =  01}  then 


0(T)  >  2  -  W(A,Xj  )  =  2  -  A1 ^  (xj  )/ [3  +  Agl  (Xl  )f.  (38) 


W(A5Xl )  has  a  single  maximum  with  respect  to  A,  denoted  by  WM (xj )  at 
A  =  3/ [2  gi(Xl)].  W„  (Xl  )  =  (4/243 )  (fj  /g,  ) ,  where  = 

2(Xl  cosh  Xl  -  sinh  Xl)/sinh  xa ,  and  WH  increases  monotonically  with 
from  its  zero  minimum  value  with  respect  to  x^^  ,  occuring  at  Xj  =0. 
For  large  X;1  ,  (fj/gj  )  «=  2 (x1  -  l) ,  and  WM  (xx  )  as  (8/243)^  ;  this  latter 
quantity  approaches  the  value  2  for  %  6l,  so  that  A  «  1049 ,  Thus, 
the  rarefaction  shock  possibility  in  the  ideal  gas  is  effectively 
ruled  out . 


B.  Single  Phase ,  Van  der  Waals 

Here  the  test  is  to  see  if  Eq.  (11 )  can  be  fulfilled.  Computations 
show  the  largest  value  of  G(V,T)  to  be  2.184  at  V/Vc  =  1.483,  ana 
T/Tc  =  .9711  (a  point  on  the  saturated  vapor  line).  Here  C° /R  ^  16.65 
will  satisfy  Eq.  (ll),  independent  of  (C°/r)':  attaining  this  condition 
requires  at  least  A  =  14  modes,  each  one  almost  completely  excited  at  a 
temperature  near  critical  for  the  substance. 


C .  Equilibrium  Coexisting  Phases 


For  this  domain,  one  searches  for  conditions  which  might  satisfy 
inequality  Eq.  (26).  Again  making  use  of  Eqs.  (34)  and  (36),  one  finds 
that  Eq.  (26)  cannot  be  fulfilled  for  A  less  than  83  vibrational  modes, 
each  of  whose  0  ^  3-5  Tc  for  the  substance.  This  restriction  for 


■pi  •v'  o  /-•+-  “1  pv  ^  . 

o  a. a.-L- 


1  purposes  nullifies  the  rarefaction  shock  probability  in  the 


equilibrium  coexistence  state  of  this  fluid. 


VIII,  comparison  with  steam  tables 

Following  the  appearance  of  the  early  analysis  by  Zeldovich  for 
the  vapor  phase  of  the  Van  der  Waals  fluid,  a  short  treatment  by 
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Novikov9  was  given  for  the  two-phase  state  of  steam.  By  using  an 
adiabatic  differential  equation  of  state  and  the  resulting  interpolating 
formula  with  (unreferenced)  steam  tables,  Novikov  concluded  that  the 
necessary  condition  for  the  rarefaction  sho  'k  could  be  satisfied  in  the 
two  phase  region  near  critical.  For  comparison,  we  show  several 
equilibrium  isentropes  for  steam  in  this  region,  using  given  tables1'1. 
This  is  done  by  noting  that  the  entropy  is  linear  in  V  on  the  isobaric 
isotherm  segment  in  coexistence;  or,  for  fixed  T  and  P, 


S(V,T)  =  S(sat.  liq.)  +  As(v-V(sat.  liq. )}/(v(sat .  liq. )  -  v(sat.  liq.)}, 


where  As  is  the  vaporization  entropy  for  this  T.  Thus,  one  can  find 
from  the  tables  a  value  of  V  for  fixed  S  and  various  T  (or  P), 

Numerical  different at ion  of  the  resulting  equilibrium  isentropes  then 

/^2  p\ 

shows  where  \gy5y  5  0  in  coexistence.  We  confirm  Novikov's  conclusion 

for  a  region  between  P/Pc  %  =  92  and  1.00.  Fig.  3  shows  an  enlarged 
view  of  the  coexistence  boundaries  for  both  steam  and  the  theoretical 
Van  der  Waals-Maxwell  fluid,  as  well  as  representative  isentropes  for 

/a2P\ 

each.  The  steam  region  where  J  <  0  is  enclosed  by  its  saturation 

locus  and  the  nebulous  hatched  boundary  indicated  at  P/Pc  ^  =9^-  This 
diffuse  lower  boundary  could  not  be  more  precisely  defined  because  of 
the  lack  of  significant  figures  in  the  steam  table  data.  All  but  one 
of  the  isentropes  for  the  theoretical  model  correspond  to  C° /R  =  3/2. 

Two  theoretical  isentropes  passing  through  the  critical  point  were 
computed,  one  for  C$/R  =  3/2  and  one  for  Cy /R  =  5/2;  this  latter  one 
has  only  a  slightly  larger  slope  than  the  former. 

Fig.  3  also  illustrates  a  common  fault  of  the  Van  der  Waals- 
Maxwell  model  in  the  critical  region:  its  scaled  saturation  locus  is 
much  less  flat  at  critical,  and  lies  well  inside  that  for  steam;  the 
same  behavior  is  noted  with  respect  to  many  other  substances16. 

Although  the  particular  steam  tables  used  here  indicate  the  necessary 
condition  for  the  rarefaction  shock  is  met  in  a  small  coexistence  region 


near  critical,  it  is  by  no  means  certain  that  su -h  a  discontinuity  would 
actually  develop.  The  critical  temperature  for  steam  is  quite  low 
{tc  (steam)  =  647°  k]  and  there  is  a  good  possibility  that  metastable 
supers aturation  of  the  liquid  phase  would  occur  during  adiabatic 
expansion.  Such  a  metastable  state  is  less  likely  at  mu  h  higher 
temperatures.  One  might  then  '■onjecture  that  if  correspondence  holds 
between  steam  and  a  fluid  with  a  much  higher  critical  temperature,  the 


If*  ■+■  -4-  -v»  <-»  ■«  i  V.  n  t-  «->  v",  /•>  --v  wi  ■?  rrVvf"  oi  v»4-  v>  i-\  -v»  rv  -f  o  n  t  -i  n  /-I  -i  o  r*i4 — i  nm  -t  t  r 

clO  OCX  ucuicc  mj.grru  buppui  o  a  iaiciauoiuii  QlblUIlLliiLii  . 


IX.  SUMMARY 


This  investigation  for  the  rarefaction  shock  condition  in  an 
equilibrium  Van  der  Waals-Maxwell  fluid  has  shown  it  to  be  practically 
impossible  in  and  across  all  phase  regions.  Comparison  with  steam 
tables,  however,  shows  that  a  necessary  condition  for  the  discontinuity 
does  exist  in  principle,  in  the  equilibrium  '-oexistence  state  near 
critical.  The  Van  der  Waals  fluid  model  used  here  is  recognized  as  not 
accurately  describing  a  real  fluid  in  the  critical  region.  Here  it 
might  prove  useful  to  perform  a  similar  investigation  using  more  modern 
concepts  about  critical  state  behavior18 .  However,  despite  promising 
progress,  these  modern  approaches  of  the  critical  region  do  not  yet 
appear  to  have  produced  a  sufficiently  simple  and  differentiable  state 
equation  as  to  enable  completion  of  such  an  investigation. 

Although  the  rarefaction  discontinuity  was  shown  to  be  impossible 
across  a  liquid-vapor  transition  boundary4,  some  solid-solid  and  liquid- 
solid  phase  transitions  do  contain  the  necessary  conditions  for  this 
discontinuity.  These  are  mentioned  in  Refs.  [4]  and  [17 1,  with  some 


experimental  evidence  for  their  existence  offered  in  Ref . 


j  • 
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APPENDIX 


The  limiting  values  of  many  quantities  making  up  the  L(vri  ,  t) 
defined  in  Eq.  (22)  for  the  coexistence  state  are  tabulated  by  Barieaue 
His  method  employed  a  power  series  expansion  of  reduced  density  along 

+  V10  hnnnHQ-mr  ■noov'  nvifi  fiol  TllP  T’PdUC'Pd  V.ipt’P 

P  =  P/Pc  ;  Y  =  T/Tc  and  ot  =  Vc/V,  and  with  j  =  1  or  3  representing 


reduced  densities  on  the  saturated  vapor  and  liquid  loci,  respectively. 
All  other  thermodynamic  quantities  for  the  coexistence  state  can  be 
expressed  in  terms  of  these  03  and  their  scaled  temperature  derivatives 
The  Maxwell  rule  gives  a  functional  relationship  between  ay  and  a3 ; 
defining  y  =  cn3  -  1  and  x  =  1  -  or  ,  the  power  series  expansion  using 


4-  Vi  Mr*. -.rr  rrs  1  1  wil  A  o 

0  lIc  r'lctAw^X-L  ruic  1^ 


■l 


a.1  x 


(Al) 


i=l 


Barieau  found  the  first  three  coefficients  to  be  ai  =1,  =  (1/5)  and 

=  (1/25).  In  order  to  find  the  limit  at  critical  of  our  functions 
B'vT),  N(t),  and  K(t)  it  was  necessary  to  find  the  fourth  order 
coefficient,  which  we  determine  to  be  34  =  (19/350).  It  is  understood 


that  within  coexistence  the  pressure,  P,  is  the  vapor  pressure  PA 


We 


further  adopt  the  notation  that  a  prime  on  a  reduced  variable  denotes 
differentiation  by  reduced  temperature,  y.  Then  our  dimensionless 
function  Z(t)  =  1  +  /B '  has  the  limiting  value  at  critical  of  8.2. 
Making  use  of  a4 ,  we  further  find  that  the  critical  limit  of  B"'  is 

-luc  iunm.iuu  u  \ij  ±S  equal  tO  /  ic  /  (.  ym  t  p  j,  a.Hu  naS  the 
limiting  value  at  critical  15.195^  (Pc/Tc*5)*  We  show  in  Fig.  4  the 
the  behavior  of  J ' (t) (Tc 3 /Pc ) ,  always  positive  in  the  range  of  interest 
The  form  of  N(t)  given  by  Eq.  (25)  does  not  permit  finding  its 
limit  simply  by  taking  the  sum  of  the  limits  of  each  term  on  the  right, 
because  these  term  limits  are  +  <=  and  -  °=;  instead,  a  tedious  expansion 
of  this  sum  of  terms  using  Eq,  (Al)  yields  at  critical  the  limiting 
value  of 


31 


~  f tvt ( m  \  1  _  fo/QUo/ 
jj.ni  c jlm v -l  y ->  —  \J /oj\x\f 


/w.  Wn  JiJihrA 


-*-0 


and  Lira  iTN(t)}  =  (3/8)r(13.444o) .  Moreover,  since  Lim  (m(t)}  =  4.5R, 
then  Lim  {k(t)/r}  =  31.860. 
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